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ABSTRACT 

Compressible flow over a laser turret causes variation in density and the 
index of refraction. As a result, a laser beam develops a phase distortion. 

Phase distortion has been calculated for both blunt and small perturbation 
turrets. For the blunt turret, the Janzen- Rayleigh technique was used to 
determine the flow field. Phase distortions of 2.2 wavelengths at 3.8 microns 
were calculated for the blunt turret. For small perturbation turrets a 
versatile analytical model was developed for a turret on a fuselage with circular 
cross section. With a two-dimensional Fourier series representation of the turret, 
any shape can be considered. Both subsonic and supersonic flows can be calculated. 
Phase distortions of 1.2 wavelengths at 3.8 microns were calculated for one turret 
at high subsonic Mach number. In addition to being of value for laser turrets, the 
methods are applicable to reconnaissance aircraft using photographic equipment and 
cruise missiles using celestial navigation. 

INTRODUCTION 

Large scale laser telescopes are being installed aboard aircraft [1, 2]. 

Airflow over the laser turret causes density variations. Since the index of 

refraction depends on density, a phase distortion is generated as the beam 
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propagates outward from the laser telescope. A phase distortion may tilt, 

defocus, or focus the beam. Higher-order aberrations, such as astigmatism 
and coma, may also occur. 

Denote the characteristic geometric scale of the laser turret by £. The 
diameter of the telescope aperture is D. The ratio D/£ tends to be much greater 
for laser installations than for other optical devices commonly installed on 
board aircraft. As an example, consider the pilot's eye. For this case, £ is 
the size of the canopy, and D is several millimeters. The pilot does not 
experience optical distortion due to the small value of D/£. 

Laser turrets tend to differ from camera installations aboard reconnaissance 
aircraft [3]. Due to the large size of D and the desire to have wide angular 
coverage, laser turrets protrude into the airstream. Camera installations usually 
have flat windows which are flush mounted in the fuselage. Furthermore, camera 
apertures may be considerably less than laser apertures. As a result of these two 
facts, degradation of photographic image quality usually is not a serious problem. 
However, in the event of image degradation, the methods of this paper could be 
applied to remove distortions. 

Cruise missiles with intercontinental ranges may use celestial navigation 
which requires precise measurement of angular location of stars. A distorted 
lens created by the ambient flow field causes an error in measurement. Knowledge 
of the external aerodynamics is important [3]. 

Optical distortion due to the external flow field can be divided into two 
categories [4]. Viscous flow phenomena fall into one category and include shear 
layers, laminar and turbulent boundary layers, and the shedding of discrete 
vortices. The other category involves the external inviscid flow field. Based 
on the preceding discussion, the reader recognizes that this paper treats the latter 
category. 
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Laser turrets can be classified as either blunt or small perturbation turrets. 
If a turret has a surface with a normal vector, n, aligned parallel or nearly 
so, to the freestream velocity vector, V , the turret is "blunt." If the 
angle between the normal vectors of the turret surface and the free- 
stream velocity is everywhere large, e.g., 60° to 90°, then the turret is 
"small-perturbation." Both types of laser turrets are discussed in this paper. 

BACKGROUND INFORMATION 
Critical Mach Number 

A useful concept is the critical Mach number, M£. At the critical Mach 
number, somewhere on the body a local Mach number is sonic. When M > M*, the 
flow becomes inherently nonlinear and shock waves appear. The transonic flow 
equations must be solved when M ot = M*. The occurrence of shock waves usually 
implies severe degradation of the laser beam quality. 

The critical Mach numbers for several body geometries are listed in 
Table I. One advantage of the solution developed later in this paper is 
the fact that the critical Mach number can be calculated. Results are shown in 
Table I for a cosine shaped turret which will be discussed in more detail in 
the section Turret on a Circular Fuselage: Subsonic Flow. 

Laser Turret Map 

Different analytical models must be developed for different laser turret 
geometries. A map can be drawn showing the various flow regions for the 
different turret geometries. Figure 1 is a laser turret map. The abscissa 
is the freestream Mach number, M^. The ordinate is the maximum slope of the 
turret. Various regions are defined by capital letters. 

The flow within region ABCD can be described with sufficient accuracy 
using the linearized small perturbation equations. The line CD represents the 
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Table I. Critical Mach Numbers for Several 

Geometric Shapes 

Body Shape 

Ratio 

,C/E 0 

Critical 

Mach 

Number 

Circular cylinder with axis normal to flow [5] 

- 

0.3985 

Sphere [5] 

- 

0.5868 

Hemisphere-cylinder with cylinder axis parallel to 
flow [6,7] 

- 

<0.7 

Cosine shaped turret on circular fuselage [3] 

0.05 

0.88 , 


0.10 

0.82 


0.15 

0.76 


0.20 

0.72 


0.25 

0.68 


0.30 

0.65 


e = amplitude or maximum height of turret 

Rg = radius of circular fuselage 

turret length/f uselage radius ratio = 1.005 


critical Mach number. The line BC is shown at a slope of 30°. The assumptions 
involved in the linearization of the potential flow equation become less and less 
valid as the maximum slope increases. The region CDJ, which resembles an inverted 
triangle, requires solution of the nonlinear, small perturbation, transonic, 
potential flow equation. 
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Figure 1. Freestream Mach Number Regions for Application of Various Analytical 
Models. 

EFGH defines the region where the solution to the blunt-body nonlinear flow 
equations is necessary. The line GH, which defines the critical Mach number, is 
shown dashed; the critical Mach number is a function of turret shape. Hence GH 
is intended to suggest qualitatively the upper bound for Mach number for 
blunt turrets. Within region EFGH the Janzen-Rayleigh [5,8,9,10] technique is 
most useful. 

Vertical lines 1, 2, 3, and 4 occur within region EFGH. These lines 
represent the upper bound of applicability of the first, second, third, and 
fourth-order solutions when the Janzen-Rayleigh technique is used. The first- 
order solution may be applied with region F11E; the second-order solution, within 
F22E; etc. The hemispherical turret considered in the following section is an 
example of a turret falling along line FG. 
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Phase Distortion, Index of Refraction and Pressure Coefficient 


The local index of refraction, n, is related to the local gas density, p , by 


n = 1 + (c^— = 1 + K ,P - 


5 

p 


SL 


( 1 ) 


where p is the freestream density at the altitude of the aircraft, and p 

00 uL 

is the density at sea level. The constant k' has a value of approximately 
2.3 x 10 ^ in the infrared. Note that K = K’Cp^pg^)* obviously < is altitude 
dependent. Assuming isentropic flow, density and pressure are functions of each 
other 


P/P m = (P/PJ Y 


Using the definition of pressure coefficient, C^, one can show that 

YM^C 


( 2 ) 


(3) 


In equations (2) and (3) y is the ratio of heat capacities, and M^ is the 
freestream Mach number. Combining equations (1) to (3) yields 


n = 1 + k 




(4) 


Since C may be positive or negative, n may be increased or decreased by the flow. 
Optical path length, L, is defined as 


b 

L ± ~ f n (s) ds (5) 

a 

where s is the distance along a ray, and points a and b are positioned on the ray. 
The index of refraction n is a function of s. The subscript i identifies a 
particular ray within the laser beam. The difference in optical path length 
between two rays i and j is 
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( 6 ) 


I 

b 

AL = Lj - L ± = J [n(s^) “ 
a 

The optical path length has the dimensions of length, e.g., meter; to make L 
nondimens ional, the wavelength of radiation, X, can be used as a reference. 

The ratio AL/A is known as the phase distortion P. Combining equations (4) and 
(6) gives a formula for P, 


1 1 



Equation (7) shows the dependence on altitude. The solutions for the various flow 
fields give values of C^. 

Zernike Polynomials 

A wavefront shape or the same thing, the phase distortion, can be expressed 
in terms of Zernike polynomials [111. Equations (8) through (17) can be found in 
the paper by Hogge and Butts [12]. 


FM = 

(dr.; 

n 

(uniform phase shift) 

(8) 

F 2 (r) = 


i 1/2 1 





1 



(9) 



»/2 | 

> (tilt) 


F 3 (r) = 


1 ,j 



(10) 

F*(r) = 

© 

1/2 / 

1 (x 

2 +y 2 ■ 

- yj , (refocus) 

(ID 


/ 6 \ 

1/2 




F s (r) = 

(r/r) 

(x 2 

- y\ 


(12) 



1/2 

i 

> (astigmatism) 


F 6 (r) = | 

(xy). 


(13) 


/ 
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( 14 ) 


^ = (^ 5 ) Ox 2 + 3y 2 - 2 « 2 k 

f, °W = (~ s)' /2 (3^ + 3/ - 2^. 


> (coma) 


(15) 

(16) 
(17) 


The Zernike polynomials are an orthonormal set of functions over an aperture of 
radius R. The phase distortion can be represented as 


10 

p = Z A j F j (r) 
j=l 


(18) 


where A^. is a constant which characterizes the phase distortion P. A_. is given by 


2tt R 


W / P(r ' 


0)F(r,6)rdrd6 


(19) 


0 0 

In equation (19) one could have used x and y as variables; however, use of 0 and r 
is more appropriate for a circular aperture. 

Since the Zernike polynomials are related to the various aberrations, 
knowledge of. the coefficients A^ indicates the magnitude of each aberration [12]. 
Further, the values of A^. are helpful for designing adaptive optics systems to 
overcome partially the aberrations. 


BLUNT LASER TURRET 

As an example of phase distortion due to compressible flow over a blunt turret, 
a hemispherical turret was selected [13,14]. Another shape amenable to analytical 
solution would be an ellipsoidal turret. 
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Janzen- Rayleigh Solution 


The geometry is shown in Figure 2. A beam of radius is propagated from a 
turret of radius R. The elevation of the beam above the y-z plane is 0 ' . 

Distance along the beam is s. Points within a plane normal to the beam are 
located by coordinates R' and a. The windward and leeward sides of the laser beam 
can be determined from the flow arrow. 



Figure 2. Geometry of Hemispherical Laser Turret Used for Phase Distortion Calculation. 


In spherical coordinates, the equations of motion are as follows: 

9uv9u_v^_ 1 Sp _ a 2 3p 

u 3r r 00 r p 3r p 3r 


3v v. 9v . uv _ 1 3p _ a 3p 

3r r 36 r ” pr 30 “ “ pr 30 


( 20 ) 

( 21 ) 


i(uf 

p 3r 


. v 3 cl 1 3 2 . 

r 30^ 2 3r ( u) 

r 


Fl'iTe !e (v sin e) = 0 


( 22 ) 
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Equations (20) and (21) are Euler's equations; equation (22) is the continuity 

equation. For flow over a sphere, the azimuth angle is not a dependent variable. 

2 2 

Multiplying equation (20) by + u/a and equation (21) by - v/a , adding, and 
combining with equation (22) yields 


1_ i_/ 2 . 1 a_ 

2 3r* r u; r sin 6 96 
r 


, . „> 1 / 2 9u , uv 9v 9v , v 9v\ 

(v sin 6) - -j (u 3 ^ + — 3 e + uv (23) 


The speed of sound, a, is the local value and changes from point to point in the 
flow. Variations in a must be accounted for. Using the energy equation along a 
streamtube, one can show that 

2 , n / 2 . 2 > 


2 2 
a 

oo 




(24) 


Inviscid flow is considered. A potential function, 4>, can be introduced for the 
velocity components 


u = U 


9 £ 

9r 


and 


v = 


_ U _9£ 


r 98 


(25) 


(26) 


Combining equations (23) to (26) yields an equation in which the only dependent 
variable is <|). 

The Rayleigh-Janzen expansion considers the potential function to be given as 

<*> = * 0 + + *& 2 + * * • (27) 

The solution is considered to be the first-order solution; the solution involving 
both and is termed the second-order solution. In this paper only the first 
and second-order solutions are discussed. For ease of writing, the following 
definition is introduced: 


* - 


( 28 ) 
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Substitution of \p into the potential equation gives 

-, 2 , . „ 2 _, w 2 


r* 0 + r-t - M‘/<i> r 4, rr + \ * r * 6 * 0r - ±3 *^ r + \ *o*b8 

\ r r r 


■) 


( 29 ) 


where all terms on the right-hand side are determined from the first-order 

2 

potential function, cj)^. The result of equating coefficients of like powers in M a 


yields two equations 


V *0 " 0 


(30) 


and 


v 2 * - mV*! - M 2 fe + \ * r Wre - *3 *ef r + \ 

\ r r r / 


(31) 


Equation (31) is Poisson's equation. 

The first-order solution can be obtained from Milne-Thompson [15], Lamb [16], 
or Karamcheti [17]. It is 


*0 ■ °( 


r cos 0 + 


R cos 0 


2r 


) 


(32) 


Equation (32) can be used to evaluate the right-hand side of equation (29). 

In July, 1916, Lord Rayleigh [10] reported the solution for the second-order 
function. The Legendre functions 


and 


= cos 0 


D 5 3 fl 3 . 

P^ = y cos 6 - — cos 0 


(33) 


(34) 


are introduced into equation (31) . The result is 
Jl 


_ u[7 36 R 6 9 R 9 \ p , L 24 R 6 3. R 9 

1 " 2[V T 7 + 2 7oj P l + [ 6 r 4 5 r 7 + 2 r 10j P 3 J 


(35) 
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Using the defining equation for the Legendre function, equation (35) can be 
solved to yield 


» - *o + - K-JKl + ™»[( ; 

+ / 3RL + 27R^_ _ 3R 1_ 3R 9 \ 1 

V 10r 2 55r 4 10r 5 176r 8 / 3 J 


R^_ 

" ' 1 3r 2 5r 5 24r 8 


h 


(36) 


Equation (36) can be inserted into equations (25) and (26) to find u and v. 

Knowing u and v, equation (24) can be evaluated. To evaluate equation (1) for n, one 
needs p/p . Since the flow is isentropic 


P_ 

P« 




1 __ 

Y-l 


(37) 


The local Mach number is given by 




(38) 


Combining equations (1) and (36) through (38) with equation (5) allows 
calculation of the optical path length. 

Geometry for Calculating Phase Distortion 

Refer once again to Figure 2. The flow properties, e.g., density and Mach 
number, are functions of r and 0. These are polar coordinates. Consequently it is 
necessary to express r and 0 in terms of 0*, s, a, and R f . The equations are 


r 2 = f (R + s)sin 0' - R 1 cos a cos 0’] 2 + R' 2 sin 2 a 

2 

+ [(R + s)cos 0* + R' cos a sin 0’] 


(39) 


and 


0 = tan 


/ 2 2 2 2 
-1 / [(R + s) sin 0* - R' cos a cos 8*] + R* sin a 

[ (R + s)cos 0’ + R' cos a sin 0'] 


(40) 
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Equations (39) and (40) allow calculation of independent variables for the flow 
in terms of a location in the laser beam as specified by s, R', a, and 0'. 

Figure 3 illustrates the distance along the beam, s, and the gap which 
exists between the surface of the turret and the s = 0 plane. The phase 
distortion is given by 



where p is the density along the axis of the beam and p is a reference density. 

U r 

The second integral in equation (41) has been termed the "gap integral." 

Graphical Presentation of Results 

A computer program has been developed for the HP9830 which calculates the 
phase distortion, P, as a function of sphere radius, R; elevation angle, 8'; radius 
within the beam, R' ; angle within the beam, a; freestream Mach number, M to ; 
wavelength. A; index of refraction constant, K ? ; ratio of heat capacities, y; 
freestream speed of sound, a^; and density ratio, p^/p^. Values used in the 
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calculations are as follows: 


beam diameter = R/2 = 0.4572 meter 
R = 0.9144 meter 
A = 3.8 x 10 ^ meter 

k' = 2 x 10~ 4 

a = 342 meters/sec 

Y = 1.4 

P «/ p SL = 1 *° 

M , R' , a, and 8’ have been varied. Calculations have been performed to the extent 
necessary to plot phase distortion maps. 

Plots of isocontours of phase shift were made for steps of 18° starting with 
0 ' = 0 and 0 ' = 90°. The calculations ignore the gap integral; evaluation of 
gap integral is somewhat arbitrary. However, to provide insight to the 
magnitude of the gap integral, the quantity was evaluated with p^_ = in 
equation (41). Results are shown in Table II. 



Table II. 

Maximum Value of Gap Integral 

Beam Elevation Angle 

0' 

Maximum Value of Gap Integral 
and Location Within Beam 


0 

0.368 

* 

18 

0.654 

0° 

36 

0.678 

0° 

54 

-0.836 

180° 

72 

-0.863 

120° 

90 

-0.866 

O 

o 


*Does not depend on a. 
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The quantity, maximum value of gap integral, is the maximum phase distortion, P, 
caused by the gap between the plane s = 0 and the surface of the turret. From 
Table II, the large value of the gap integral indicates that the gap cannot be 
ignored . 

The plots of isocontours of phase shift are shown in Figures 4 to 9. To avoid 
awkward decimal values, the phase distortion has been multiplied by 100. Hence 
250 from the graphs is 2.5 wavelengths phase distortion; a value of 25 from the 
curves would be a quarter wavelength. Recall the beam axis is used as the 
reference for the phase shift, P. Reference to equation (41) shows that P is 
positive when p along the ray in question exceeds p^. A positive value of P 
means the wavefront lags behind the front at the beam axis. 

The critical Mach number 0.587 was chosen for the calculations. The outer 
edge of each plot is, of course, the outer edge of the beam. The ratio of 
turret diameter to beam diameter is 4.0. The significance of the maps will be 
discussed in the following subsection. 



Figure 4. Contours of Constant Figure 5. Contours of Constant 

Phase Shift for 8' = 0°. Phase Shift for 0' = 18°. 


E 









An important question is the rate of decay of the integrand in equation (41) . 


Calculations were made of 


3P 

3 (s/R) 

which is the phase distortion per unit distance along the beam. Results of the 
calculation are shown in Figure 10. The phase distortion, P, is the area 
enclosed by one of the curves. Most of the phase distortion occurs within a 
distance of one turret radius. 



RATIO OF DISPLACEMENT ALONG BEAM 
TO TURRET RADIUS, */ R 

Figure 10. Phase Distortion per Unit (s/R) as a Function of (s/R) for Rays 

at a = 0°, 60°, 120°, and 180° Within the Beam. Elevation angle was 
0’ = 18°. 


The windward side of the beam is at 0 = 0°. When 0’ = 18 and 0=0, the 
air is compressed and P is positive. When 0 = 180°, i.e., on the leeward side 
of the beam, the air is expanded and P is negative. 

Interpretation of Results 

The contours in Figure 4 are a series of concentric circles with a negative 
value of phase shift. Hence, the beam is being focused. The focal length, F, can 
be estimated from 


F = 



(42) 
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Using R^ = R/4 = 0.2286 meter, A = 3.8 micron, and P = 1.6, the value of F is 
8.6 km. 

For 0' = 54° in Figure 7, the contours are almost straight lines. The 
algebraic signs of phase shift indicate the beam is being tilted in a direction 
opposite to the relative wind, i.e., the beam is leaning into the wind. The tilt 
angle is given by 

rJ P 

(tilt angle) = A -^r (43) 


Inserting values from Figure 7, the tilt angle is found to be 26.5 microradians. 

This tilt angle should be compared with the pointing accuracy of the pointer- 
tracker which forms the turret. Fortunately tilt can be corrected easily by 
adaptive optics. 

For 0’ = 90°, the flow over the turret causes defocusing of the beam; 
see Figure 9. Using equation (42), the focal length is F = - 6.3 km. The flow gives 
effectively a negative lens. 

Scaling of Phase Distortion 

Using equation (41) one can demonstrate the scaling relationships for altitude 
(density), turret size, and laser wavelength. The equation is 


p i X 2 p i R i K i 


(44) 


where subscripts 1 and 2 refer to two different hemispherical turrets at the 
same M^, The phase distortion becomes more severe as wavelength decreases. 
Throughout the infrared region, the value of k' is nearly constant; however, as 
one approaches the visible, <' tends to increase with decreasing wavelength. Of 
course, for constant diffraction, the ratio R/A will tend to be constant. Since p 
follows an exponential decay with altitude, the phase distortion falls rapidly 
with increasing altitude. 
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TURRET ON A CIRCULAR CROSS SECTION FUSELAGE: SUBSONIC FLOW 


Using the solution for the wavy wall on a cylinder [18] and Fourier 
analysis, a turret of any shape can be described mathematically, and the flow 
field can be obtained. A question about this model concerns the effect of fuselage 
ends. 

Fuselage End Effects 

For the analytical model to apply, the fuselage shown in Figure 11 must extend 
to infinity in both directions. Obviously real aircraft have finite length. What 
is the influence of aircraft fuselage ends? 




Figure 11. A Small Perturbation Turret on a Circular Cross Section Fuselage. 


To gain insight to this question, the pressure distribution on the surface 
of several axisymmetric bodies has been calculated assuming potential flow. The 
fineness ratio of the fuselage, i.e., the ratio of length to diameter, was varied. 
The potential equation for axisymmetric flow is 


aV 


8z 


,2 


+ L. ML) = o 

+ 3r ' ' r s-'' u 


3r' 


(45) 
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where z' is distance in the flow direction and r* is radial distance using 
cylindrical coordinates. Equation (45) assumes incompressible flow although the 
calculated flow may be accurate to = 0.3, depending on body bluntness. 

A body surface is generated by a source-sink combination separated by 
a distance £. Lengths are made nondimensional using £. The potential 
function is nondimensionalized by the following 

4> - 4>7 u/ (46) 


where is the frees tream velocity. 

There are two other lengths in the problem in addition to £. First, there 
is the body diameter. Second, there is the distance from the source to the 
upstream stagnation point which is designated as z*. Once again z* is a 
nondimensional length using £ as a reference. Incidentally, the upstream and 
downstream stagnation points are saddle point singularities; this is a fact 
to be remembered when the equations are integrated to obtain the body surface. 

The velocity components are u = u’/U^ in the z-direction, and v = v'/U 
in the r-direction. The equations for u. and v are as follows: 

_ -1 . Z* 2 (z* - l) 2 z _ Z* 2 (z* - 1) 2 (Z - 1) , 

u l + 3/2 3/ 2 

(1 - 2 z*)(z 2 + r 2 ) (1 - 2 z*)[(z - l ) 2 + r 2 ] 

and 


v = 


2 2 
z* (z* - 1 ) r 


z* (z* ~ 1 ) r 


3/2 3/2 

(1 - 2 z*)(z 2 + r 2 ) (1 - 2 z*)[(z - l ) 2 + r 2 ] 


(48) 


The flow is in the direction of the positive z axis with the source at z = 0 
and the sink at z = 1. Consequently, the upstream stagnation point occurs 
at a location where z* is a negative number. When z* approaches negative 
infinity, the body approaches a sphere. When z* approaches zero, the length 
of the fuselage, which is L = £(1 + 2z*) , approaches infinity. Any value of z* 
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satisfies equation (47) when u = 0 (a stagnation point), z = z*, and r = 0. 
However different values of z* change the length to diameter ratio of the body. 

To find the body surface, the differential equation relating velocity 
components and the slope of the body surface is solved numerically; the equation 
to be solved is 


dr ii 
dz v 


(49) 


Integration starts at z = - |z*| and r = 0. Since the rear stagnation point is a 
saddle point singularity, one cannot integrate to z = + |z*| and r = 0. The flow 
is known to be symmetric relative to z = 0.5; this fact was used to obtain the body 
for z > 0.5. 


The pressure coefficient in terms of nondimensional velocities is 

P " P °° 2 2 
= l-u - V 


Cp P U '/2 


(50) 


At the front and rear stagnation points, has a value of unity. When is zero, 
the local static pressure exactly equals the static pressure at infinity, and 
the local flow velocity is U^. When is negative, the local flow velocity 
exceeds U . 

OO 

Results of sample calculations are shown graphically in Figures 12 to 17 . 

In Figure 12 the length to diameter ratio, L/D, is 1.004, which is nearly a sphere. 

The curves from A to B and from E to F are the pressure coefficient along the 

stagnation streamlines. The curve BCDE is on the surface of the body. For a 

sphere, the pressure coefficient is 

C = 1 - f sin 2 0 (51) 

P 4 

When 0 = 90°, C = - 5/4. When 0 = 41.8°, C =0, which corresponds to point C in 
P p 

Figure 12. 

For all figures, the sharp positive peaks are at = + 1.0. Note that for 

L/D = 1.004, the C curve for the body surface is concave. For L/D = 2, the C 
P P 
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Figure 16. C for L/D = 8. 
P 


Figure 17. C for L/D = 10. 
P 


curve is convex at z = 0.5; this is point G in Figure 13. Although not shown, for 
L/D = 1.58, the curve dips negative but remains flat for 0.1 < z < 0.9. In 
Figure 13, the x's are the source and sink. In each figure there are two vertical 
lines passing through the source and sink; the vertical lines are located at 
z = 0 and z — 1. 

As L/D Increases, it is apparent from the figure sequence that at z = 0.5 

becomes less negative and approaches zero. At L/D equal to 6, C is less than 

P 

3 per cent of the positive peak for the stagnation point. A laser turret installed 
at z = 0 on a fuselage with L/D = 6 would experience freestream velocity within 
1.5 per cent. Using equation (50) and the binomial expansion, one can show that 

(52) 

Application of equation (52) to the case of L/D = 6 gives u = 1.015. 


u 


C 

1 2 . 

2 
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The Boeing 707-320 has a fineness ratio of 10; and the French Caravelle, 10.4. 
The stretch DC-8 has L/D considerably more than 10. The Boeing 747SP, which is a 
standard 747 with a section of fuselage removed, has L/D near 6. Another wide body 
jet, the DC-10, has L/D near 8.6. The laser turret on the Airborne Laser 
Laboratory, ALL, is on a fuselage with L/D of 10, and the turret is located at 
nondimens ional z = 0.3 approximately. Looking at Figure 17, one notes C^ is 
nearly zero at z - 0.3. Consequently, the turret should be experiencing near 
freestream velocity. (The preceding statement ignores interference effects of 
other components, e.g., wing.) 

A potential flow solution provides insight to the effect of fuselage ends on 
the airflow likely to be present near a laser turret. For L/D = 6, a turret 
near z = 0.5 will have essentially freestream velocity. For L/D = 10, freestream 
conditions prevail in the interval 0.3 < z < 0.7. 

Formulation of Subsonic Flow over ji Laser Turret 

For small perturbation flow, the potential equation as derived in 
Liepmann and Roshko [18] is 

+ $ 2 <p + <j> + - <p + ~ <p QQ = 0 (53) 

r 

The (+) sign applies to subsonic flow; and the (-) sign, to supersonic flow, 
the definition of (3 is 

8 2 = |l - M 2 | (54) 

The function <j> is the perturbation potential yielding only the changes from the 
freestream velocity. Equation (53) is given in cylindrical coordinates which are 
appropriate for the cylindrical fuselage. A solution to equation (53) can be 
obtained using the separation of variables. For the flow field external to the 
fuselage 

<j>(r,8,x) = A cos n8 sin qx K (gqr) (55) 
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where A is a constant, n is an integer, and q is to be determined. See Fuhs and 
Fuhs [2] for development of the solution in detail. 

The turret shape is defined by 

R(x,6) = R q + eX(x)T(6) (56) 

where X(x) and T(0) are functions. The height of the turret is e as shown in 
Figure 11. As an example, a cosine shaped turret will be used. For a cosine 
shaped turret 

X(x) = -|(1 + cos (57) 

and 

T(0) = ~(1 + cos f0) (58) 

Equations (58) and (59) apply for |x| _< £ and 1 0 [ _< 2n/f. Outside this region 
R = Rq. The quantity 1/f is the fraction of the 2TT-circumference occupied 
by the turret, and £ is the length of the turret. For a cosine shaped turret, 
the appropriate Fourier series are 


and 


X(x) 


e£ ejL ^ sin pitt£/L 1 mTrx 

2L L ^ mir£/L , , 0/TN 2 COS L 


(59) 


T (0) 



z 

n=l 


sin mr/f 
nw/f 


— cos n0 

1 - (n/f) Z 


(60) 


In equation (59) the spacing between periodic turrets on an infinitely long 
fuselage is L. The turret generated by equations (56) , (59) , and (60) is shown 
in Figure 11. 
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At the body surface for inviscid flow, there can be no flow through the wall. 
Mathematically this is expressed by 

n • V = 0 (61) 

where n is a unit vector normal to an element of the surface and V is the vector 
representing the local fluid velocity. If one can express the shape of a body by an 
equation of the form f(r,8,x) = 0, then 


n = $f/|$f| (62) 

Using the potential function for determining V, the boundary condition at the body 
surface is 

= 0 (63) 

where $ is the full potential related to the perturbation potential by 

$ = ip + xU (64) 


For small perturbation, equation (63) reduces to 

I + i = ° 

Equation (65) is derived in Appendix A of reference [3]. 

Solutions for Subsonic Flow 

In view of equations (56) , (59) , and (60) , the potential function is assumed to 
have the form 


<J>(r,0,x) = I I <j> (r,0,x) 

r n m T nm v ’ * 

Using appropriate nondimensional variables, <j> is found to be 


2A cos n0 sin K (B^iTr/L) 
I / q \ _ nm L n 

W r ’°’ x; ~ “ erK n+1 (enm/L) + K n _ 1 (3mTT/L) ] 


K n is a Bessel function discussed in the book by Hildebrand [19]. 


A is given by 
nm 6 


( 66 ) 


(67) 

The constant 
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. _ e& P 1 sin(mTT&/L)~ }{ ~ 1 sin(mr/f )~ 1 f( . R . 

“ " " fL Ll - (WL ) 2 W/L) JLl - (n/f) 2 Wf) J (68> 

Combining equations (66) to (68) , one has the complete perturbation potential 

equation. The perturbation velocities are calculated using V<f> ; these velocities 

are shown explicitly in reference [3]. With velocities in nondimens ional form, 

the pressure coefficient for axisymmetric flow is 

C = - 2u - v 2 (69) 

P 

Equation (69) can be evaluated, and the results inserted in equation (7) to 
calculate the phase distortion P. 

Before discussing typical results for phase distortion, calculation of 
critical Mach will be presented. Equation (3) relates p/p^ to the pressure 
coefficient, C^. Assume the flow is isentropic; this assumption permits use 
of the following relation between static pressure where the local Mach number is 
unity, p*, and the static pressure at infinite distance from the turret, p^: 



(70) 


Combining equations (3) and (70) establishes an equation for the critical 
pressure coefficient 


* 

C 

P 




(71) 


When C given by equation (69) equals C*, the local velocity is sonic. Equations 
P P 

(69) and (71) were used to calculate the critical Mach numbers given in Table I. 
Typical Results 

As an example, the phase distortion was calculated for a laser beam leaving the 
turret of Figure 11. The beam was pointed at an angle 90° to the axis of the 
fuselage and was symmetrically located relative to the turret. The beam radius was 
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A/2. The pressure coefficient and optical path length were calculated for a 
ray starting at x = - A/2 and for the reference ray at x = 0. The pressure coefficient 
along the two rays is shown in Figure 18. The values used in the calculation were 
as follows: 

M ra = 0.62 R q = 1 meter L = 5 

£ = 0.35 A = 3.8 microns 0=0 

p /p =1.0 A = 1.005 K ' = 0.00023 

b L 

The quantity 3P/3(s'/Rq) was evaluated for the two rays specified above. 

The results are shown in Figure 19, which is a graph of increment of phase 
distortion as a function of radial distance s'. About one-half of the phase 
distortion occurs within a distance along the beam of approximately 0.35 Rq. 

Expressed in terms of turret height, E, one-half of the phase distortion is 
generated within approximately one turret height. 

The integrated phase distortion, which is the area under the curve, is 1.17. 

In terms of wavelengths, the integrated distortion is 1.17 A. 



s'/R 0 s'/R 0 

Figure 18. Pressure Coefficient Figure 19. Increment of Phase Distortion 

along Two Rays Selected along Beam, 

for Phase Distortion 
Calculation. 
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TURRET ON A CIRCULAR CROSS SECTION FUSELAGE: SUPERSONIC FLOW 

The potential flow equation for supersonic flow is equation (53) with use 
of the (-) sign. Once again the turret shape is illustrated by Figure 11. A more 
general formulation for turret geometry is developed in this section. 

Turret Geometry 

The function X(x) of equation (56) is a polynomial 


K 


^ k 
a k x 


X(x) = 1 + X 
k=l 

Likewise the function T(0) of equation (56) is another polynomial 

P 

T(0) = 1 + ^ ^0 J 

j=2,4. 


(72) 


(73) 


Note that only even powers of 0 are used since the turret is symmetric in the 
0-direction. As before, £ is the height of the turret at x = 0 = 0. Each term in 
the polynomial is expanded separately as a Fourier series [20]. 


* k -Z 

m=0 


mTTX , , . IMTX 

a . cos — — + b , srn — 
mk L mk L 


Hence 


X(x) = X 


mirx . , . mTTx 1 

V cos ~C * b „k sln IT 


k=l \m=0 


In a similar manner 


(74) 


(75) 


0^ = a . cos n0 
^ nj 


n=0 


(76) 
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Consequently 


T ( 0 ) = £ f Z a nj C0S n6 ] 

j=2,4, \n=0 


(77) 


The Fourier coefficients are as follows: 

ik 


& r , n k , . mirfl, kL , 

a = — [ (- 1) +1] sxn — b , , 

mk B L mTT m,k-l 


. - ^ r t i \ k+1 1 1 nnr£ . kL 

mk mTT L mTT m,k-l 


(78) 

(79) 


Equations (78) and (79) are recursion relations. To start the sequence one 


needs a„, and b„, which are 
Ok Ok 


= [(- l) k + lH k+1 
a 0k 2L(k + 1) 


(80) 


and 



(81) 


Also both a .. and b n are zero. Beginning with k = 0, the Fourier 

nij “J. 

coefficients for any order polynomial can be obtained using equations (78) to 

(81). If k is an even number, then b ^ is zero. If k is an odd number, then 

a , is zero, 
mk 

Equations (72) through (81) define a turret of length 2£ which is periodic 
every x = 2L. In the circumferential direction, the turret is positioned 
between - 0^ <_ 0 _< 0 , where 0^ defines the boundary of the turret. To evaluate 
equation (77), one needs a n j * Equations (78) to (81) can be used to determine a^. 
by making the following substitutions: 

m -*• n L tt 

k -»■ i i -*■ 0 

J t 

Note that j has only even values. 
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not all coefficients may be independent. Due to equations (56), (72), and (73), the 
turret height at x = 8 = 0 is always e. To have a continuous surface for the 
fuselage, X(x) = 0 for |x| = £; also, T(0) = 0 for | 0 | = 9^. In addition, 
conditions may be imposed such that dX/dx = 0 at |x| = SL, etc. 

For each boundary condition imposed, one coefficient of the polynomial is 
eliminated as an independent variable. The coefficients which are not used to 
satisfy boundary conditions can be used as design variables. In a companion 
paper [21], the design variables are selected so as to optimize phase distortion. 
Solution t£ Potential Equation for Supersonic Flow 

Once again separation of variables is used to solve equation (53) . The 
solution has the form 

<j>(r,0,x) = R(r)0(0)X(x) (82) 

where 

R(r) = AJ (Bqr) + BY (Bqr) (83) 

n n 

0(0) = C sin n0 + D cos n6 (84) 

X(x) = E sin qx + F cos qx (85) 

Both J and Y are Bessel functions discussed by Hildebrand [19]. For symmetry with 
n n 

respect to 0 = 0, C is zero. Since both J and Y are finite for large r, neither 

n n 

A nor B can be set equal to zero. In contrast to the subsonic flow solution, the 
supersonic flow solution has one additional constant to be determined. Matching 
waves in the far field permits evaluation of the additional constant [22]. For 
a simple cosine wavy wall on a cylinder, the potential equation has the form as 
derived by Fuhs [23] 

<j)(r,x) = A[(sin ax - cos ax)Y^(aBr) + (sin ax + cos ax)J^(aBr)] (86) 
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The constant A is related to the amplitude of the cosine wavy wall. Equation (86) 
implies an eigenvalue given by 


Y 1 (a8R Q ) = - J 1 (a3R 0 ) 


(87) 


The wavenumber of the wavy wall is a. 

The requirement for no flow through the turret boundary is given by 
equation (65) . For the turret of Figure 11, the boundary condition becomes 


v _ Sr\ _ 1^ 9(j>\ 

U 3x/__ r " U drj 

r_R 0 r R 0 


( 88 ) 


The surface of the turret is given by 

R(0,x) = R + V Ta cos + B sin ^ X " l cos n0 (89) 

0 Z-i L nm L nm L J 

n m 

Based on equation (86) , the potential function is assumed to have the form 


4>(r,0,x) = U «, x + X Z^nm 


(90) 


n m 


where 


! , n n , „f/ . nmx mirx \„ /8m7r \ 

a (r,0,x) = A cos n0 (sin — cos — ; — lY I— — r) 

nm nm |_\ L L / n\ L / 

(91) 

+ (sin — + cos 9J 

Equation (87) establishes an eigenvalue, which is the root of equation (87). 

For a Fourier series, the wavenumber a is determined by the series. The flight Mach 
number determines 8. Hence is not arbitrary. As a result one must introduce 
an eigenvalue cylinder which has the same wavenumber as the cylinder required by 
the Fourier series; however, the eigenvalue cylinder has a different radius. 
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Further, the eigenvalue cylinder is shifted in phase. The equation for a 
cosine shaped eigenvalue cylinder is 

R (0,x) = R + e cos n0 cos — (92) 
nm mn mn L 

The boundary condition for the eigenvalue cylinder is 



nm 


(93) 


The phase shift, x , is given by 

x' = x + x (94) 

nm 

Carefully note that R (0,x) and d) are formulated with x' whereas R(0,x) is in 

nm nm 

terms of x. 

Consider first the cosine wavy wall defined by equation (92). Use of the 
boundary condition of equation (93) gives the eigenvalue equation 


3mirR 


nm 


f Y „+i <c »> - 


n-t-i. m n m n m 


(95) 


where 


BirmR 


C = 
m 


nm 


(96) 


There is only one value of R which satisfies equation (95); hence equation (95) 

nm 

provides a method to calculate R in equation (92) . 

nm 

From the boundary condition of equation (93) , one finds 

X nj c J Y n+1 ( c J + J n +i< c J] - n[Y (C ) + J(C >]} 
nm l m n+1 m n+1 m n m n m • 


nm 


mTTR 


(97) 


nm 


Using the boundary condition for the turret as given by equation (88), the phase 

shift, x , can be determined 
nm 
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tan 


mux D - E 

nm _ run nm 

L D + E 

nm nm 


(98) 


where D and E are defined as 
nm nm 


firmr BlTmR BtHDR 

D = - -E-— y ( ~) + Y ( ^ 

nm L n+l ' 1 L } V L 


and 


(99) 


)3mTrR 


BrnnK 


nm 


Ml j ( ZI0, n_ /rro , 

L Vl ( L )+ R V 


( 100 ) 


Continuing with the boundary condition for the turret, define 


and 


a = D + E 
nm nm nm 


b = D - E 
nm nm nm 


( 101 ) 


( 102 ) 


The amplitude of <f> term, which is A , is related to A of equation (89) bv 
Iim nm nm N J 


m77A 


nm 


nm 9 0 1/2 

L[a Z + ] 

nm nm 


(103) 


The value of is determined by turret shape. Combining equations (91), (94), 
(98) , and (99) through (103) , one obtains 


4> (r,6,x) = - 

nm 


(mir/L)A cos n9 , 

nm > % mTTx „ ^Bmirr, 

(a + b ) sin — — Y (-=- ) 

nm nm L n v L ' 


2 2 

a + b 
nm nm 


, _ , s mTTx „ ^Bnnrr* 

- (a - b ) cos — — Y (-— ) 

nm nm L n L 


+ Ka - b > sin SS + (a + b ) c03 JH]j (SHE, 
nm nm L ' nm nm L n L 


(104) 


Following a similar procedure for the sine shaped wavy wall, one finds 
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<(> (r,0,x) 

ran’ = - 


(mTr/L)B cos n0 
ran 

2 

a + b 
ran ran 


✓ it \ mux „ ,8nnrr N 

** (a__ + b ) cos — - Y (— — ) 
ran ran L n L 


. / , . . mux „ .fhnirr,. 

+ (a - b ) sin — Y (— ) 

ran ran L n L 


r / , . mirx , . \ j nnrx-j T , Bmirr. 

- [(a - b ) cos — (a + b ) sin — ]J (-^ ) 

run nm L ran nm L n L 


The combined potential equation is 


(105) 


<{)(r,0,x) = U^x + ^ <>nm + (106) 

n m 

Equation (106) can be used to determine the velocity components. With the 
velocity components, equation (69) can be evaluated for the pressure coefficient, 
0^. In turn, the pressure coefficient can be inserted into equation (7) to 
determine the phase distortion, P. 

Results of ^ Sample Calculation 

The flow over a cosine shaped turret was determined. The turret was 
defined by equations which approximate a cosine 

X(x) = 1.0 - 0 . 5 (x/&) 2 + 0.0625(x/£) 4 (107) 

and 

T(0) = 1.0 - 1.82(0/0 t ) 2 + O.O832(0/0 t ) 4 (108) 

The Mach number was = 2.0, e = 0.2, H = 2Rg, and 0 - 60°. 

Figure 20 shows the phase distortion map for the case being considered here. 
The distortion is almost pure tilt. 
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Figure 20. Map of Phase Distortion. Azimuth 0°; Elevation 45°; Mirror 
Radius = R^/20. 

CONCLUSIONS 

Laser turrets may be classified as blunt or small perturbation depending 
on maximum slope of the turret surface. Analytical models were developed 
for both categories of turrets. 

As an example of a blunt laser turret, the flow over a hemispherical 
turret was solved using Janzen- Rayleigh expansion technique. Terms to the 
second order were obtained and were used to calculate the phase distortion. 

At the critical Mach number for a hemispherical turret, the phase distortion 
was several wavelengths. The phase distortion varied from a positive focus for 
0° elevation to combined focus and tilt at intermediate elevations and to 
negative focus for 90° elevation. The distortion can be expressed in terms of 
Zernike polynomials. 
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A versatile analytical model was developed for a laser turret on an 
infinitely long circular cylinder. The laser turret is described by a two-dimensional 
Fourier series using the solution for the flow over a wavy wall on a cylinder. The 
turret can have any shape subject to the constraint that the turret is a small 
perturbation in the flow. Both subsonic and supersonic solutions were obtained. 

End effects were examined. For a fuselage with length to diameter ratio of 
6.0, the laser turret mounted at the midpoint of the fuselage would experience 
freestream velocity with 1.5 per cent. For an aircraft with L/D of 10, such as 
ALL, any location between z/D =0.3 to 2/D = 0.7 would have freestream conditions. 

The analytical model permits calculation of critical Mach number for the 
turret. The critical Mach number is important because of the change from 
linearized subsonic flow solution to nonlinear transonic flow. 

For the subsonic case, the maximum observed distortion was a fraction of 
wavelength. Similar results were found at Mach 2. 

The supersonic flow solution has eigenvalues which require introduction 
of an eigenvalue cylinder. A cylinder of arbitrary radius cannot yield a 
satisfactory solution. A cylinder slightly smaller can satisfy the boundary 
conditions and becomes the eigenvalue cylinder. For large radial distances from 
the cylinder, the solution approaches that of a planar wavy wall. 

In addition to useful calculations for laser turrets, the analytical models 
are useful for reconnaissance aircraft using cameras and cruise missiles using 
celestial navigation. 
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